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^ Abstract 

> 

This paper introduces QuSAnn vl.2 and Multiplexor Expander vl.2, two Java ap- 
plications available for free. (Source code included in the distribution.) QuSAnn 
is a "code generator" for quantum simulated annealing: after the user inputs some 
parameters, it outputs a quantum circuit for performing simulated annealing on a 
quantum computer. The quantum circuit implements the algorithm of Wocjan et 
al. (arXiv:0804.4259), which improves on the original algorithm of Somma et al. 
>. (arXiv:0712.1008). The quantum circuit generated by QuSAnn includes some quan- 

tum multiplexors. The application Multiplexor Expander allows the user to replace 
each of those multiplexors by a sequence of more elementary gates such as multiply 
controlled NOTs and qubit rotations. 



1 Introduction 

For an explanation of the mathematical notation used in this paper, see some of my 
previous papers; for instance, Ref.pQ Section 2. 

We say a unitary operator acting an array of qubits has been compiled if it 
has been expressed as a SEO (Sequence of Elementary Operations, like CNOTs and 
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single-qubit operations). SEO's are often represented as quantum circuits. 

There exist software, "general quantum compilers" (like Qubiter[2]), for com- 
piling arbitrary unitary operators (operators that have no a priori known structure). 
There also exists software, "special purpose quantum compilers" (like each of the 7 
applications in QuanSuite[3l IH E]), for compiling unitary operators that have a very 
definite, special structure which is known a priori. 

This paper introduces QuSAnn vl.2 and Multiplexor Expander vl.2, two Java 
applications available for free. (Source code included in the distribution.) QuSAnn 
is a "code generator" for quantum simulated annealing: after the user inputs some 
parameters, it outputs a quantum circuit for performing simulated annealing on a 
quantum computer. QuSAnn is not really a quantum compiler (neither general nor 
special) because, although it generates a quantum circuit like the quantum compilers 
do, it doesn't start with an explicitly stated unitary matrix as input. Multiplexor Ex- 
pander is not a quantum compiler either for the same reason. Multiplexor Expander 
can be more aptly described as a "code translator": it takes a SEO and replaces it 
by a different but equivalent SEO. 

In Ref. [6], Somma et al. proposed an algorithm for quantum simulated an- 
nealing that requires, for any e > 0, order 1/Vo elementary operations to find, with 
probability greater than 1 — e, the minimum of a function. Here 5 is the distance be- 
tween the two largest eigenvalue magnitudes of the transition probability matrix for 
the Metropolis Markov chain. The algorithm of Somma et al. outperforms the clas- 
sical simulated annealing algorithm, which requires order 1/5 elementary operations 
to do the same thing. 

Subsequently, Wocjan et al. in Ref. [7] improved on the algorithm of Somma 
et al. (See also Refs.[SJ[H], where Wocjan et al. discuss related issues). 

Both the Somma et al. and the Wocjan et al. algorithms use Szegedy quantum 
walk operators[10] and phase estimation. But the Wocjan et al. algorithm uses a 
Grover fixed point search [11] instead of the quantum Zeno effect. 

The quantum circuit generated by QuSAnn implements the algorithm of Woc- 
jan et al. given in Ref. [7J. The circuit includes some quantum multiplexors. (See 
Ref. [lj for a review of quantum multiplexors.) The application Multiplexor Expander 
allows the user to replace each of those multiplexors by a sequence of more elementary 
gates such as multiply controlled NOTs and qubit rotations. Multiplexor Expander 
gives the user the option of expanding the multiplexors in two different ways: either 
as an Exact SEO (see Refs. [2J[T]), or an Oracular Approximation (see Ref. [12"]). 

A nice feature of the source code for QuSAnn, Multiplexor Expander and all 
the applets in QuanSuite, is that they all share a common Java class library (named 
QLib). 

2 QuSAnn 

The QuSAnn applet makes the following 3 assumptions: 
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1. For its annealing schedule (see Appendix |A| for definition) /3q, fli, fa, ■ ■ ■ Pt f , if 
assumes fa = and fa + i — /3j = A/? > for all j = 0, 1, . . . , tf — 1. 

2. For its energy function (see Appendix [A] for definition), it assumes E(x) = 
(x — ^f-) 2 , where the state space (i.e., the set of x values) of the minimization 
problem is {0, 1, 2, . . . , N s - 1}. 

3. For its neighborhood function (see Appendix |A| for definition), it assumes neig(x, y) = 
1 if \x — y\ < 1 and neig(x, y) = otherwise. 

These 3 assumptions were made in order to make the applet simple. They can 
be easily changed (i.e, one can use a more complicated annealing schedule, energy 
function and neighborhood function) by making trivial alterations to the source code 
of QuSAnn. 

2.1 The Control Panel 

Figjl] shows the Control Panel for QuSAnn. This is the main and only window of 
the application. This window is open if and only if the application is running. 
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Ver. 1.2 
Inputs 



QuSAnn 



...... 

...... 



QuSA nn 



File Prefix 

Number Of State Bits 
Number Of Probe Bits (for each PE step) 
Number Of Phase Estimation (PE) Steps 
Crover Depth 
Upper Bound on Number of Neighbors 
Number of Betas fz~ 
Delta Beta Per Unit Time 



nip 



Output Highlights 
Number of Qubits 
Number of Elem. Ops. 
Message 



( Write Files ) 



Figure 1: Control Panel of QuSAnn 



The Control Panel allows you to enter the following inputs: 

File Prefix: Prefix to the 3 output files that are written when you press the Write 
Files button. For example, if you insert test in this text field, the following 3 
files will be written: 

• test_qsann_log.txt 

• test_qsann_eng.txt 

• test_qsann_pic.txt 



Some examples of these output files are given in Section 2.2 below. 



Number of State Bits: The parameter Ng = 1,2,3... defined in Appendix|X} The 

state space of the minimization problem is {0, 1, 2, ... , 2 Nb — 1}. The transition 
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probability matrix M of the Metropolis Markov chain is a 2 Nb dimensional 
matrix, and the Szegedy quantum walk operator W(M) is a 2 2Nb dimensional 
matrix. 

Number of Probe Bits (for each PE step): The parameter a = 1,2,3,... de- 
fined in Appendix [Dj See Figfl8| 



Number of Phase Estimation (PE) Steps: The parameter c = 1,2,3,... de- 
fined in Appendix [Dj See FigJT8] 



Grover Depth: The parameter df = 1,2,3,... defined in Appendix D df is the 
final level of recursion to which one wishes to carry out the fixed point Grover 
search. 

Upper Bound on Number of Neighbors: The parameter upBdNeig £ M >0 de- 
fined in Appendix [A| 

Number of Betas: The parameter tf+ 1 = 2, 3, . . . defined above to be the number 
of betas in the annealing schedule. 

Delta Beta Per Unit Time: The parameter A/3 £ 1R >0 defined above to be the 
difference between adjacent betas of the annealing schedule. 

The Control Panel displays the following outputs: 

Number of Qubits: The total number of qubits used by the circuit, equal to 2Nb + 
ac in the notation of Appendix |D} 

Number of Elementary Operations: The number of elementary operations in 
the output quantum circuit. If there are no LOOPs, this is the number of lines in 



the English File (see Sec. 2.2.2), which equals the number of lines in the Picture 
File (see Sec. 2.2.3). For a LOOP (which is not nested inside a larger LOOP), 
the "LOOP k REPS : N" and "NEXT k" lines are not counted, whereas the lines 
between "LOOP k REPS : iV" and "NEXT k" are counted N times. Multiplexors 
expressed as a single line are counted as a single elementary operation (unless, 
of course, they are inside a LOOP, in which case they are used repeatedly). 

Message: A message appears in this text field if you press Write Files with a bad 
input. The message tries to explain the mistake in the input. 



2.2 Output Files 

Figs. [2j [3]and[4j were all generated in a single run of QuSAnn (by pressing the Write 
Files button just once). They are examples of what we call the Log File, English 
File, and Picture File, respectively, of QuSAnn. Next we explain the contents of 
each of these output files. 
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f \ /"\ ^ test_sann_log.txt 

Inputs : 

File Prefix = test 
Number of State Bits = Z 
Number Of Probe Bits (for e 
Number Of Phase Estimation 
Grover Depth = Z 
Upper Bound on Number of Ne 
Number of Betas = Z 
Delta Beta Per Unit Time = 

Outputs : 

test_sann_eng .txt 
test_sann_pic .txt 
test_sann_log . txt 
Number of Qubits = S 
Number of El em. Ops. = 1096 
A 



Figure 2: Log File generated by QuSAnn 

2.2.1 Log File 

Fig|2] is an example a Log File. The Log File records all the information found in the 
Control Panel. 

2.2.2 English File 

Fig|3] is an example of an English File. The English File completely specifies the out- 
put SEO. It does so "in English", thus its name. Each line represents one elementary 
operation, and time increases as we move downwards. 

In general, an English File obeys the following rules: 

• Time grows as we move down the file. 

• Each row corresponds to one elementary operation. Each row starts with 4 
letters that indicate the type of elementary operation. 

• For a one-bit operation acting on a "target bit" a, the target bit a is given after 
the word AT. 

• If the one-bit operation is controlled, then the controls are indicated after the 
word IF. T and F stand for true and false, respectively, al stands for a control 
n(a) at bit a. aF stands for a control n(a) at bit a. 

• "LOOP k REPS : iV" and "NEXT k" mark the beginning and end of iV iterations, 
k labels the loop, k also equals the line-count number (first line is 0) of the line 
"LOOP k REPS:iV" in the English file. 
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BY -0.0 


mp_y 
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IF 
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3T 




BY -0.0 


mp_y 


AT 
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IF 
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4(0 


3T 




BY 


-60.00000000 
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LS0.0 


i.T 
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BY 
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6(0 


4(Z 


ZT 




BY -0.0 














.r .n. 









Figure 3: English File generated by QuSann in the same run as the Log File of FigfJ} 
Bottom of file is not visible. Right hand side of file is not visible. 



• SWAP a (3 stands for the swap (exchange) operator E(a,/3) that swaps bits a 
and (3. 

• PHAS 9 degs stands for a phase factor e l9de3S ^o. 

• POPH 9 de9S stands for the one-bit gate e iP ° 0dess ns . P1PH 9 de9S stands for the 
one-bit gate e tPl6 e9S iio. Target bit follows the word AT. 

• SIGX, SIGY, SIGZ, HAD2 stand for the Pauli matrices ax, Cy, o"z and the one-bit 
Hadamard matrix H, respectively. Target bit follows the word AT. 

• ROTX, ROTY, ROTZ, ROTN stand for one-bit rotations with rotation axes in the 
directions: x, y, z, and an arbitrary direction n, respectively. Rotation angles 
(in degrees) follow the words ROTX, ROTY, ROTZ, ROTN. Target bit follows the 
word AT. 

• MP_Y stands for a multiplexor which performs a one-bit rotation of a target bit 
about the y axis. Target bit follows the word AT. Rotation angles (in degrees) 
follow the word BY. Multiplexor controls are specified by a(k, where integer a 
is the bit position and integer k is the control's name. 

Here is a list of examples showing how to translate the mathematical notation 
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Figure 4: Picture File generated by QuSAnn in the same run as the Log File of Fig|2} 
Bottom of file is not visible. 



used in Ref. [T] into the English File language [j] 



1 This same table appeared before in Ref. [3] except that we have added a new row of boxes at 
the end of the table for multiplexors expressed as a single line in the English and Picture files. 
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Mathematical language 


English File language 


Loop named 5 with 2 repetitions 


LOOP 5 REPS: 2 


Next iteration of loop named 5 


NEXT 5 


£(l,0f (3)n(2) 


SWAP 1 IF 3F 2T 


e i42.7 I f 5 n(3)n(2) 


PHAS 42.7 IF 3F 2T 


e i42.7 T | 5 n(3)n(2) 


POPH 42.7 AT 3 IF 2T 


e i42.7 T |jn(3)n(2) 


P1PH 42.7 AT 3 IF 2T 


(Tx(1) n(3)n(2) 


SIGX AT 1 IF 3F 2T 


ay(1) n(3)n(2) 


SIGY AT 1 IF 3F 2T 


az (l)«(3)n(2) 


SIGZ AT 1 IF 3F 2T 


H H\n{3)n(2) 


HAD2 AT 1 IF 3F 2T 


f e l T f 5 23.7<rx(l))n(3)r i (2) 


ROTX 23.7 AT 1 IF 3F 2T 


( e | T f 5 23.7<7y(l))n(3)n(2) 


ROTY 23.7 AT 1 IF 3F 2T 


( e 3lfe23.7«Tz(l))n(3)fi(2) 


ROTZ 23.7 AT 1 IF 3F 2T 


(etiii[ 30,7 ^( 1 ) +40,Ty ( 1 ) +11<TZ ( 1 )])"( 3 ) n ( 2 ) 


ROTN 30.0 40.0 11.0 AT 1 IF 3F 2T 


where < 


9 blbo o- Y (3)P blbo (2,l)in{0) 

' 9 00 = 30.0(^) 
e i = 10.5(^) 
flio = 11.0(^) 

k ^ii = 83.1(j^) 


MP_Y AT 3 IF 2(1 1(0 OT BY 30.0 10.5 11.0 83.1 



2.2.3 ASCII Picture File 

Fig|4] is an example of a Picture File. The Picture File partially specifies the output 
SEO. It gives an ASCII picture of the quantum circuit. Each line represents one 
elementary operation, and time increases as we move downwards. There is a one-to- 
one onto correspondence between the rows of the English and Picture Files. 
In general, a Picture File obeys the following rules: 

• Time grows as we move down the file. 

• Each row corresponds to one elementary operation. Columns 1, 5, 9, 13, . . . rep- 
resent qubits (or, qubit positions). We define the rightmost qubit as 0. The 
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qubit immediately to the left of the rightmost qubit is 1, etc. For a one-bit 
operator acting on a "target bit" a, one places a symbol of the operator at bit 
position a. 

• I represents a "qubit wordline" connecting the same qubit at two consecutive 
times. 

• -represents a wire connecting different qubits at the same time. 

• + represents both I and -. 

• If the one-bit operation is controlled, then the controls are indicated as follows. 
@ at bit position a stands for a control n(a). at bit position a stands for a 
control n(a). 

• "LOOP k REPS : iV" and "NEXT k" mark the beginning and end of N iterations, 
k labels the loop, k also equals the line-count number (first line is 0) of the line 
"LOOP k REPS : N" in the Picture File. 

• The swap (exchange) operator E(a,f3) is represented by putting arrow heads < 
and > at bit positions a and f3. 

• A phase factor e l ° for 6 G R is represented by placing Ph at any bit position 
which does not already hold a control. 

• The one-bit gate e lPo ^ e for 6 G R is represented by putting OP at bit position 
a. 

• The one-bit gate e tPl ^ e for 6 G R is represented by putting @P at bit position 
a. 

• One-bit operations ax(cx), 0y(a), cr^(a) and H(a) are represented by placing 
the letters X,Y,Z, H, respectively, at bit position a. 

• One-bit rotations acting on bit a, in the x, y, z, n directions, are represented by 
placing Rx,Ry,Rz, R, respectively, at bit position a. 

• A multiplexor that rotates a bit r about the y axis is represented by placing 
Ry at bit position r. A multiplexor control at bit position a and named by the 
integer k is represented by placing (k at bit position a. 

Here is a list of examples showing how to translate the mathematical notation 
used in Ref. p] into the Picture File language]^] 

2 This same table appeared before in Ref. [5] except that we have added a new row of boxes at 
the end of the table for multiplexors expressed as a single line in the English and Picture files. 
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Mathematical language 


Picture File language 


Loop named 5 with 2 repetitions 


LOOP 5 REPS: 2 


Next iteration of loop named 5 


NEXT 5 


£( 1;0 )n(3)n(2) 


— @ — < — > 


e i42.7 I f 5 n(3)n(2) 


@ +— Ph 


e i42.7 T |jn(3)n(2) 


OP — @ I 




e i42.7 T f f; n(3)n(2) 


OP — @ I 




(Tx(1) n(3)n(2) 


o— o— X 




ay (!)n(3)«(2) 


@ Y 




az (l)W(3)n(2) 


o— o— z 




_ff(l)«(3)n(2) 


@ H 




( e llfo 23 - 7CT x(l)W(3)n(2) 


@ Rx 




^1^23.7(7^(1)^(3X2) 


@ Ry 




^ e | I f 5 23.7a z (l)\n(3)r l (2) 


@ Rz 




( e l ii 5 [30^(l)+40 t ry(l) + llaz(l)])n(3)n(2) 


@ R 




|^g* Z/6i,60 

where < 


9 6l( , ov(3)n i 6 (2,l)]n(0) 

' ^oo = 30.0(^) 
e i = 10.5(^) 

10 = 11.0^) 

011 = 83.1(jfo) 


I Ry— (1— (0— Q 



3 Multiplexor Expander 

QuSAnn outputs a quantum circuit which includes multiplexor operations. Multi- 
plexor Expander can read the output files of QuSAnn and write new files in which 
each multiplexor is replaced by a sequence of more elementary operations such as 
multiply controlled NOTs and single qubit rotations. Multiplexor Expander gives 
the user the option of expanding the multiplexors in two different ways: either as an 
Exact SEO (see Ref. [21 1]), or an Oracular Approximation (see Ref. [T2"]). 
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3.1 The Control Panel 



Fig|5] shows the Control Panel for Multiplexor Expander. This is the main and 
only window of the application. This window is open if and only if the application is 



running. 



O ^ rA 



Multiplexor 



Multiplexor Expander 



Ver. 1.2 
-Inputs — 

Prefix for Input Files 
Prefix for Output Files 
Compilation Mode 
Bit Precision 



Output Highlights 
Number of El em. Ups. 
Message 



Exact SEO 



r 



U3 



(Write Output Files j 



Figure 5: Control Panel of Multiplexor Expander 



The Control Panel allows you to enter the following inputs: 

Prefix for Input Files: Prefix to the 2 input files that are read when you press the 
Write Output Files button. For example, if you insert testln_qexp in this 
text field, the following 2 files will be read: 

• testln_qexp_eng.txt 

• testln_qexp_pic.txt 



Some examples of these input files are given in Section 3J2 below. These 2 files 
are usually English File and Picture File outputted by QuSAnn. If they aren't, 
they must be formatted in the same way as the 2 files outputted by QuSAnn 
or else Multiplexor Expander will fail. 

Prefix for Output Files: Prefix to the 3 output files that are written when you 
press the Write Output Files button. For example, if you insert testOut_qexp 
in this text field, the following 3 files will be written: 
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testOut_qexp_log . txt 
testOut_qexp_eng . txt 
testOut_qexp_pic . txt 



Some examples of these output files are given in Section 3.3 below. 



Compilation Mode: The compilation mode, either as Exact SEO (see Ref. [21 H]), 
or Oracular Approximation (see Ref. [T2]). 

Bit Precision: The number of significant fractional bits in the oracular approxima- 
tion (see Ref. [12]). This parameter is ignored if the compilation mode is Exact 
SEO. 

The Control Panel displays the following outputs: 
Number of Elementary Operations: Same as in QuSAnn Control Panel. 
Message: Same as in QuSAnn Control Panel. 

3.2 Input Files 

Figs. [6] and [7] are examples of the 2 input files for Multiplexor Expander, what we 
call the Input English File and Input Picture File, respectively, of Multiplexor 
Expander. These examples are not really output files of QuSAnn, but they are 
formatted in the same way as English and Picture files of QuSAnn. 

© O r> _ testln_qexp_eng.txt 



MP_Y AT 4 IF 3(Z 2(1 1(0 0T BY -0.0 -154.843Z5773)8 



Figure 6: Input English File for Multiplexor Expander. Right hand side of file is not 
visible. 



eee 


j>_ testln_qexp_pic.txt 


Ry-(Z-(1- 


(0— @ 





Figure 7: Input Picture File for Multiplexor Expander. 
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3.3 Output Files 



Figs. [8j [9] and 10 were all generated in a single run of Multiplexor Expander (by 
pressing the Write Output Files button just once), with 

• the input files of Figs. [6] and [7| and 

• the compilation mode set to Exact SEO . 



Figs. [8j [9] and 10 are examples of what we call the Output Log File, Output 
English File, and Output Picture File, respectively, of Multiplexor Expander. 
The notation of these files is the same as that for the Log, English and Picture files 
for QuSAnn (see Section 2.2). 



Figs. 11, 12 and 13 are similar to Figs. [8j [9] and 10 above, except that they 
were all generated in a single run of Multiplexor Expander, with 

• the input files of Figs. [6] and [7| and 

• the compilation mode set to Oracular Approximation. 



©O© L~l te5tOutExact_qexp_log.txt 


Inputs : 




Prefix for Input Files = testln_qexp 




Prefix for Output Files = testOutExact_qexp 




Compilation Mode = EXACT_SEO 




testln_qexp_eng .txt 




testln_qexp_pic .txt 




Outputs : 




testOutExact_qexp_eng . txt 




testGutExact_qexp_pic .txt 




testOutExact_qexp_log . txt 




Number of Elem. Ops. - 16 






/a 



Figure 8: Output Log File generated by Multiplexor Expander with Exact SEO as 
compilation mode. 
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testOutExact_qexp_eng.txt 



ROTY 


351.S9459Z7SZ7Z756 


AT 


4 


IF 


0T 


SIGX 


AT 4 IF 


IT 


m 






ROTY 


-59.39459Z78Z7Z7635 


AT 


■1 


I r 


ST 


SIGX 


AT 4 IF 


ZT 


m 






ROTY 


-S1.S9459Z7SZ7Z7&4 


AT 


■i 


I" 


ST 


SIGX 


AT 4 IF 


IT 


err 






ROTY 


-Z10.&05407Z17Z7ZZ4 


AT 


■i 


IF 


ST 


SIGX 


AT 4 IF 


ST 


9T 






ROTY 


S1.S9459Z7SZ7Z7&4 


AT 


■1 


IF 


ST 


SIGX 


AT 4 IF 


IT 


ST 






ROTY 


-59.39459Z7SZ7Z7635 


AT 


■1 


I r 


ST 


5IGX 


AT 4 IF 


ZT 


ST 






ROTY 


-S1.S9459Z7SZ7Z7&4 


AT 


■1 


I r 


ST 


SIGX 


AT 4 IF 


IT 


ST 






ROTY 


59.39459Z7SZ7Z766 


AT 


■1 


I" 


ST 


SIGX 


AT 4 IF 


3T 


ST 







Figure 9: Output English File generated by Multiplexor Expander in the same run 
as the Log File of FigJH] 



8 6 8 li testOLitExact_q ex p_pic.txt 




Figure 10: Output Picture File generated by Multiplexor Expander in the same run 
as the Log File of FigM 
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6 O 6 testGutGrac_qexp_log.txt 


Inputs : 




Prefix for Input Files = testln_qexp 




Prefix for Output Files = testOutOrac_qexp 




Compilation Mode = ORACULAR_APPROX 




Bit Precision = 3 




testln_qexp_eng . txt 




testln_qexp_pic . txt 




Outputs : 




testGutOrac_qexp_eng .txt 




testOutOrac_qexp_pic .txt 




testOutOrac_qexp_log .txt 




Number of Elem. Ops. = 13 









Figure 11: Output Log File generated by Multiplexor Expander with Oracular Ap- 
proximation as compilation mode. 







teatOutOrac 


_qexp_eng.txt 






ROTY 


360.0 


AT 


4 


IF 


ST 


0T 






SIGX 


AT 


6 


IF 




2- 


IT 


0T 




SIGX 


AT 


e 


IF 


3" 


2T 




0T 




SIGX 


AT 


6 


I r 


3^ 


IT 


IT 


0T 




SIGX 


AT 


e 


I r 


3T 


IT 


IT 


0T 




ROTY 


1S0.0 


AT 


A 


IF 


fcT 


ST 






SIGX 


AT 


e 


IF 


3" 




IT 


0T 




SIGX 


AT 


t 


I r 


3^ 


ZT 


_ r 


0T 




SIGX 


AT 


6 


I r 


3^ 


ZT 


IT 


0T 




SIGX 


AT 


6 


I r 


3T 


ZT 


IT 


0T 




SIGX 


AT 


7 


IF 


3T 


ZT 


: r 


0T 




ROTY 


90.0 


AT 


4 


IF 


7T 


m 






SIGX 


AT 


7 


I r 


3T 


ZT 




0T 























Figure 12: Output English File generated by Multiplexor Expander in the same run 



as the Log File of Fig, 11 
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© © [sj testOutOrac_qexp_pic.txt 



1 


1 §— Ry- 


-+ — + — +— 




1 


X— +— +— 


_a— o— g>— 




1 


X— +— +— 


-0— g— 0— 


-e 


1 


X— +— +— 


-0— 




1 


X— +— +— 


-@— i— e~ 




1 


@- --+--- Ry- 


-+ — + — + — 


,X 


1 


X— +— +— 


_0— o— g— 




1 


X— +— +— 


_0_— ft— 0— 


@ 


1 


X— +— +— 


-o— e— @— 




1 


X— +— +— 


-@---g---g— 




X 


— + — + — + — 


_gi g o — 




?: 


--+---+- --Ry- 


— h — + — + — 


a- 


X 


--+ — + — +— 







,4 



Figure 13: Output Picture File generated by Multiplexor Expander in the same run 



as the Log File of Fig, 11 



A Appendix: Classical Simulated Annealing 

The goal of simulated annealing, both classical and quantum, is to find the minimum 
of a (bounded below) function (i.e. to solve a minimization problem). The function 
E : S x — > R-° to be minimized will be called the energy function (assumed non- 
negative without loss of generality). It's domain S x will be called the state space. 
We will assume that Sx = {0, 1,2, .. . , Ns — 1}, were Ns = 2 Nb is the number of 
states and Nb is the number of bits. Besides the state space and energy function, it 
is also convenient to introduce a neighborhood function neig : Sx x Sx — » Bool, 
defined by neig(x,y) = 6{x and y are neighbors). 

To solve this minimization problem, classical simulated annealing uses a Markov 

chain. 

A Markov chain is a Bayesian network x —>■ x_ x — > • • ■ —>■ x t , wherein 
all random variables G S^ have the same range of values: Sx. = Sx for all 
j, and every node except the first one has the same transition probability matrix: 
Px +1 \x- (y\ x ) = M(y\x) for j — 0, 1, . . . , tf — 1 and x, y £ S^. A stationary state 7rQ 
of the Markov chain with transition probability matrix M is a probability distribution 
on Sx which is also an eigenstate of M with unit eigenvalue, 

M( V \x)tt(x) = 7r(y) (1) 
for all y £ S x - We say that a probability distribution 7r(x) is a detailed balance of 
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M if 

M(y\x)7r(x) = M(x\y)7r(y) (2) 

for all x, y £ S 1 ^. Clearly, if 7r() is a detailed balance of M, it is also a stationary state 
of it. 

Classical simulated annealing uses a special Markov chain due to Metropolis. 
The Metropolis transition probability matrix Mp for a given minimization problem 
and inverse temperature (3 > 0, is defined as follows: 



(x ^ y) 



neig(x,y) 



f3[E(y)-E(x)} 



} 



Mp(y\x) = { \ '"7'^ ^ , (3) 

where upBdNeig is some real number greater or equal to max y ^2 X neig(x, y). Thus, 



upBdNeig is an upper bound on the number of neighbors. Fig|T4] tries to explain the 
logic behind Eq.([3]). A "system" prefers going downhill to going uphill, but is willing 
to visit a neighbor living uphill occasionally. 



i-: 



neig(x, y) e -p(Efy)-E(x» < 1 
jpDdNeig — upBdNeig 



Figure 14: Transition probabilities for Metropolis Algorithm. 



One can show that the following probability distribution (called a normalized 
Boltzmann factor) is a detailed balance and therefore a stationary distribution of 

p -f3E(x) 

*p{*) = • (4) 

Zp is called the partition function. It is defined so as to make 7r() a probability 
distribution: 

x 

In classical simulated annealing, we consider a product Mp t . . . Mp 1 Mp of 
transition probability matrices. The monotonically increasing (or at least non-decreasing) 
sequence of non-negative real numbers (3q, (3\ . . . , (3 tf is called the annealing sched- 
ule. 
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B Appendix: Q-Embedding of Probability Matrix 



In this Appendix, we will review the concept of q-embedding of a probability matrix, 
as used in Ref . [13] and later in Ref . [H] . 

Given a conditional probability P(y\x) where x G S ¥ and y G S y , we will call 
a probability matrix the matrix P with entries P(y\x) with rows labeled by the y 
and columns by the x. Any unitary matrix U with matrix elements (y\(x\U\y)\x) = 
A(y, x\y,x), where x,x G Sx and y,y G S y , that satisfies 

J2\My^\y = o,x)\ 2 = P(y\x), (6) 

X 

for all a;, y, will be called a q-embedding (quantum-embedding) of the prob- 
ability matrix P. AQ acts like a probability amplitude. The index y that we set 
to a fixed value (call it zero) is called a source and the index x that we sum over is 
called a sink. This nomenclature is similar to the one used by Fredkin and Toffoli 
when they showed how any boolean function / : Bool 171 — > Bool n can be embedded in 
a reversible function. See Ref. [13] for more details and references. Note that Eq.(|6]) 
is satisfied if we set 

A(y,x\y = 0,x) = 6^P(y\x) . (7) 

A q-embedding of a probability matrix is of course not unique. Next we 
will give one possible q-embedding for any square probability matrix acting on Nb 
bits. The q-embedding that we will give is very convenient because it is expressed 
succinctly as a product of quantum multiplexors. (See Ref. [1] for a review of quantum 
multiplexors). 

We begin by pointing out some trivial algebraic results that will be used below. 
Note that for any 9 G R, 



Co —Se 
So Cf> 



Cg, Se are shorthand for cos 9 and sin^, respectively. Thus, if b G Bool, we get from 
the first column of this matrix that: 

f C e if 6 = 

(b\e-^\0) = C b 9 S b 9 ={ . (9) 

I 5^ if 6 = 1 

(We are using the notation of Ref. [1] for b where b G Bool, namely = 1 and 1 = 0.) 
A compact expression for the entries of both columns can be obtained as follows. For 
a, b G Bool, 

(b\e- iaY& \a) = (b\e- taY V a |0) = (b © a\e~ i{ - 1)a(TYe \0} . (10) 
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The last matrix element in Eq.(lO) is given by Eq.([9]). 

Let q e M> NsxNs be the probability matrix for which we desire a q-embedding. 
First consider Nb = 1. We want the following constraint to be satisfied: 



A(bi,b Q \ax = 0,a ) = ^ v / #iKj 
where bi, b , at, a e Bool. Define a unitary matrix U by 



11' 



U 



where the angles 9\ ao are defined by 




Lj a On 

V\a Q eVo 



9l = »il 



a ■ 



Using P b \b') = 5%\b) for b, b' e .Boo/, and Eqj9| it follows that 



:i2i 



(13) 



5»o| — C> — |oo) 
>i|-CH0> 



^ a o °(&i|e" iCTy(1)e|a °|0) 



bQ V\a V\a 



Hence U is a q-embedding of q. 

Now consider Nb = 2. We want 



A(b 3 ,b 2 ,b u b \a 3 = 0,a 2 = 0,ai,ao) = 8%8%y/q(b3, hWi, oo) , 
where all dj and bj are in Bool. Define a unitary matrix U by 



(14a) 

(14b) 
(14c) 



(15) 



U 




- i J2b 2 ,a 1 ,a cr Y( 3 ) e b 2 \a 1 a Pb 2 a 1 a ('2,i-fi) -iJ2 ai ,a °Y ( 2 ) | a X a Q Pay a a (1 ,0) 



where the angles 0& 2 [ oiao and 9\ aiao are defined by 



Q2b 3 g2& 3 Qb3b2\aiao 



(16) 
(17) 



and 
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Cgf* Sg b2 — q. b 2 \a 1 a ■ (18) 

(We are using the notation of Ref . [TJ where a dot at the position of an index means 
that the index has been summed over; e.g., q.b = q a b)- It follows that 




C K° (h\e- iaY (3)e "2 i«i«o |0) (b 2 \e~ icrY {2)e ^ -o |0> (1 9a) 



°bi °b 



r^b 3 qb 3 

O/i On 

b 2\ a l a b 2\ a l a 



On On 



(19b) 
(19c) 



Thus, as in the Nb = 1 case, U is a q-embedding of q. 

It's clear how to generalize this construction so as to get a q-embedding of a 
probability matrix q G M NsxNs for any positive integer Nb- 



C Appendix: Szegedy Quantum 
Walk Operator W 

In this appendix, we will review the definition and some useful properties of the 
Szegedy quantum walk operator W (first defined by Szegedy in Ref. \±0\, first used 
for simulated annealing by Somma et al. in Ref. j6], first implemented in terms of 
multiplexors here). 



C.l Symmetric Matrix M sym 

For any Markov chain with transition probability M, define matrix A (the entry-wise 
square root of M) by 

K{y\x) = VM(y\x) , (20) 
and the matrix M sym (a symmetric version of M) by 

M sym {y\x) = A(x\y)A(y\x) , (21) 

for all i,i/6 S^.(Note that unlike M(y\x), M sym (y\x) is not a probability function in 
y, its first argument.) 

Define the quantum states 

i(*r> = x>r(aom ( 22 ) 

X 
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for T] = |, 1. (Note that only the rj = | state is normalized in the sense of quantum 
mechanics.) 

Claim 1 

M\n) = |tt) , (23) 

and 

M sym |V^) = |v^> • (24) 
t4/so, M and M sym have the same eigenvalues. 

proof: 

Taking the square root of both sides of the detailed balance statement Eq.Q, 

we get 



K{y\x)^Ax) = A{x\y)y^y) . (25) 

Therefore, 



M, 



sym 



(y\x) = A(x\y) r -^ A{x\y)^y) = —±=M(x\y)y/n(y). (26a) 



'7r(x) V 7r ( a; ) 

Hence, 

M(y\x)7r(x) = ^ M(x\y)n(y) = n(y) , (27) 

X X 

and 

VM s ^XyN)v / ^ = ^^^^(^b)\/^Mv / ^R= V^M- (28) 

x x ^/< X ) 

Order the elements of the finite set S x in some preferred way. Use this preferred 
order to represent M and M sym as matrices. Define a diagonal matrix D whose 
diagonal entries are the numbers ir(x) for each x G Sx, with the x ordered in the 
preferred order: 

D = diag[(ic(x)\ x ] . (29) 

Since 



i , , „ i 



M sym =D~*MD3, (30) 
it follows that 
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det(M - A) = det(M sym - A) (31) 

for any A G R. 
QED 

Let the eigenvalue^] of M sym (and also of M) be m , m 1; . . . m^-i G K with 
m = 1 > |mi| > |m 2 | . . . > |mjv s -i|. Define \rrij) to be the corresponding eigenvec- 
tors of M sym (but not necessarily of M). Thus 

M sym \mj) = mj\mj) , (32) 

for j = 0, 1, ... , Ns — 1. In particular, |m ) = | v^")- 

For each j, define ipj G [0, |] and ^ G {0, 7r} so that = e tr,j cos^j. (Thus, 
cos</9j > and e ir,j = ±1). Note that m = 1 so y?o = 0. The M eigenvalue gap 5 is 

2 

defined as 5 = 1 — |mi|. 5 ~ ^ when is small. 
C.2 Q-Embeddings C/ and U 

Next we will consider two "dual" disjoint sets of distinct qubits with Nb qubits in 
each set. Let these two set be labeled a = (ai, a 2 , . . . , QLn b ) an d (3 = (Pi, Pi, . . . , Pn b )- 
Any operator Q acting on the tensor product of a state \x)s and a state where 
x, y G BooI Nb , can be represented, depending on taste, either in quantum circuit 
notation or Dirac notation, by 

o ° = ^fe/)^)* . (33) 



Let | denote the operator that swaps all bits aj and (3j for j = 1, 2, . . . , Nb- 
Let U be any unitary matrix satisfying for any x G BooI Nb , 

''''' or ^|0)|x) = (A|ar))|ar) . (34) 



U 



|0> H^r - l x > 



Let 



A(y',y\0,x) = ivl ^r {X) = (y'\(y\U\Q)\x) . (35) 

(y'K-i- |o> 



for x,0,y,y' G BooI Nb . Then, by virtue of Eq.(34) 



3 Thcrc must be a single eigenvalue 1 and all others must have a magnitude strictly smaller than 
one because of the Frobenius-Perron Theorem. The eigenvalues must be real (but they can be 
negative) because M sym is a Hcrmitian matrix. 
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A(y',y\0,x) = 5(y,x)A(y'\x) = 8(y,x)^M(y'\x) 

Thus, U is a q-embedding of the probability matrix M. 
If we define U by 



(36) 



U =\ U I , (37) 

then we can immediately write the following equations, which are dual to equations 
we wrote previously for U : 



U 



|0> 



A 



\x) 



or C/|x)|0) = |x)A|x) , 



x 



A(y',y\x,0) 



(y\- 



|0> 



(y1(y|£>k>|0) , 



i(y',y|x,0) = 5(y',x)A(y\x) = 6(y',x)y/M(y\x) 
Next define the unitary operator U by 

U = U ] U . 

Clearly, 

Note that [/ | is Hermitian and its square equals one: 

(u \y =\u ] = u\ , (u if = i . 

Matrix U has the following highly desirable property: 
Claim 2 For any j, k e {0, 1, . . . , N s - 1}, 



(0| rr |m fc ) 

[/ = mi" 



Pi 



|0> 



proof: 



(0| m \m k ) 
( mj \ |0) 





' (y|A T " 




\x)~ 


mj\y) 


(y\ 




A\x) 



E 



(x|m fc ) 



(38) 

(39) 
(40) 

(41) 

(42) 

(43) 



(44) 



(45a) 

(45b) 
(45c) 
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QED 



C.3 Definition of W 



Now define the projection operator tx by (expressed below in 3 alternative but equiv- 
alent notations) 



7T 



|0)(0| 



lf N " ®P® Nb =P (a) . 



(46) 



Define the projection operator n dual to n, by 



7T =t 7T 1 = 



|0)(0| 



(47) 



Define a reflection operator (—1)* (expressed below in several equivalent notations) 
and its dual reflection operator (— If as follows: 



-1)* = 1 - 2tt = 4 



.ifo NB = (_i)n^^) ? 



(48) 



(-If =1 (-If I • (49) 
Finally (gasp!), we are ready to define the Szegedy quantum walk operator 



W corresponding to the transition probability matrix M, by 

W(M) = U(-1)*U*(-1)* . 

C.4 Eigenvalues of W 

To find the eigenvalues of W , we will use the following identities. 
Claim 3 

7r|mj0) = \rrijO) , 



(50) 



(51a) 



n(U 1)\rrij0) = m,j\mj0) , 

for allje{0,l,...,N s -l}. 
proof: 

From the definition of n, we see that 

. |0> |0) 



7T 



K> \ m i) 



(51b) 



(52) 
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Also, 



10} 



|0)(0| u \mj) 



QED 



\rrij) k \m k )(m k \ |0) 



An immediate consequence of Claim [3] is that 



10} 

\ m j) 



(53) 



{mjO\U I \m k 0) = (mjO\TcU J \m k 0) = rrijSj 

for j,A;e{0,l,...,JVs-l}. 

Note that since m = 1, Eq.(54) implies that 



(54) 



|m 0} = C/ | |m 0} . (55) 

Another consequence of Claim [3] is that |moO} is a stationary state of W. 
Indeed, one has 



W\m 0) = U(-lfW{-lf\m o 0) (56a) 

= C/|(l-27r)|C/t(-l)|m Q 0) (56b) 

= (l-2mo^I)(-l)|moO) (56c) 

= (l-2)(-l)|mo0> (56d) 

= |m 0} . (56e) 



Let 



Vl usy = span{\mfi),Ul\mfi)} (57) 

for j G {0,1,..., Ng — 1}. (By "span" we mean all linear combinations of these 
vectors with complex coefficients.) 

Claim 4 WV 3 buS y = VLy for all j G {0, 1, . . . , N s - 1}. For j = 0, let 

l^o} = \m 0) . (58) 
{|"0o}} is an orthonormal basis for V® usy and W\iI)q) = \ipo). For j ^ 0, let 

fe> = nt? (e-^U I |m,0> - |m,0» . (59) 
is an orthonormal basis for V 3 busy and W\i/j±j) = e ±l2tpi \i)±j). 

proof: 

Using the identities of Claim [3j one finds after some algebra that 
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and 



for all j. 



W\mjO} = (-l)|mj-0> + (2m j )U J \rrijQ) 



W(U DlrrijO) = {-2 mj )\ mj 0) + (-1 + Am))U \ \mfi) 



(60a) 



(60b) 



According to Eqs.(60), Vl usy is invariant under the action of W for each j. By 



virtue of Eq.(|54|), V busy is 1-dim for j = and 2-dim if j ^ 0. We've already proven 



that |?tiq0) is a stationary state of W. 



e 2j >A e^iUl^O) 




h°H e y> 

Figure 15: Definition of \e\j) and \e2j). 



Now consider the case j ^ 0. Both U{— iyU^ and (— l) n are reflections in the 
planar subspace V^, , and reflections are a special type of rotation about the axis 
normal to this plane, so their product is also a rotation about this axis. The vectors 
\mj0), and U | |m_j-0)} are independent but not orthogonal. However, we can express 
them in terms of orthogonal vectors (see Fig, 15) as follows: 



and 



\mj0) = \e Xj ) , 



e I \rrij0) = cos(^)|e w ) + sin(^-)l e 2j) • 



(61a) 



(61b) 



W 



(62) 



In the \eij), \e 2 j) basis, we find after substituting rrij = e ir,J cos(^) into Eqs.(60) that 

cos(2tfj) sin(2<^j) 
— sin (2^-) cos(2(fj) 
The eigenvalues of this matrix are e w , with corresponding eigenvectors: 

1 



(63) 



These eigenvectors satisfy 



(64) 
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By expressing \e\j) and \ezj) in Eq.(63) in the original basis, we get Eq.(59). 
QED 

Define the following vector spaces: 

V = span{\x) (g) \y) : x,y G S^} , 



(65) 



Va = span{\x) ® |0) : x G S^} 



(66) 



V 



(67) 



and 



V 6us3/ = V A + V B ■ (6* 
V can be expressed as a direct sum of Vbusy and its orthogonal complement V^ sy : 



V = V< 



busy 



V, 



busy ' 



From Claim [4, it follows that Vbusy is a direct sum of the subspaces V\ 



•3 

busy' 



N s -1 

Vbusy = Vbusy 
j=0 



(69) 



(70) 



Recall that matrices M and M sym are iY^ dimensional whereas W is iVj dimensional. 



Since the size of is N s , dim(V) = N§. From Eq.(70) and Claim [IJ dimiVbusy) = 
2Ns — 1. Furthermore, (IV'i) : J = 0, ±1, ±2, . . . , ±(Ns — 1)} is an orthonormal basis 

for Vbusy 

At this point we've explained the action of W on Vbusy, but we haven't said 
anything about the action of W on V^ sy . Next we show that W acts simply as the 
identity on V^ sy . (This is what one would expect since the vectors in V^ gy are parallel 
to the axis of the W rotation.) Recall that if S and T are subspaces of a vector space 
V, then (S + T) L = S ± n T L . Therefore, 



From the definitions of Va and Vb, it's easy to see that 

Vj_ = span{\x) ® \y) : x G S^, and y e S^ - {0}} 

and 



(71) 



(72) 



vi = I (vj 



(73) 
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Claim 5 



for all |0) G Vj 
proof: Let 

\0) e VI 



busy ' 



busy 



w\ 



Vj n Vi. Hence 10) G Vj and 



(74) 



[/ | |0) for some 



C/(— l)*C/ t (— 1)*| 



c/ : (-1)* i c/t(-i)°i0> 

U l (1 - 2tt) J tf+tf | |0> 
C/ t (1 - 2tt)|0> 

10)- 



(75a) 
(75b) 
(75c) 
(75d) 



QED 

C.5 Multiplexor Implementation of W 

Consider the case N B = 2. Using e m = — 1 and |0)(0| —n, one gets 



w = u{-ifu ] {-\y 



u 



(76) 



Now we need to find a SEO for the U in Eq.(76). Using Eq.(41) to express U in 



terms of U, and using the method given in Appendix |B| for implementing U in terms 
of multiplexors, we get 




(77) 



In Eq.(77), a box with a dagger in it represents the Hermitian conjugate of the box 



without a dagger and acting earlier on the same qubit. It's clear how to generalize 
this construction of W to any number iV^ of bits. 



Figs{16] and 17 show English and Picture files, written using the format of 
QuSAnn and Multiplexor Expander, for a Szegedy quantum walk operator W, for a 
case with = 2. 
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PHA5 


180.0 


AT 


z 


IF 


IF 


0F 


MP_Y 


AT 


S 


I r 


3(1 


Z(0 


-91.094809947173; 


mp_y 


AT 


: 


I" 


'K- 


Z(0 


0(Z -76.570Z- 


mp_y 


AT 


j 


I r 


2(2 


1(1 


0(0 -7S.570Z- 


mp_y 


AT 


z 


I" 


1(1 


0(0 


-91.094809947173: 


PHAS 


180.0 


i.T 


1 


IF 


3F 


ZF 


mp_y 


AT 


Z 


I" 


1(1 


0(0 


91. 094809947173 8' 


mp_y 


AT 


3 


IF 


2(2 


1(1 


0(0 76.570Z4 


mp_y 


AT 


1 


IF 


3(1 


Z(0 


0(Z 76.570Z4 


mp_y 


AT 


Z 


I" 


3(1 


Z(0 


91. 094809947173 8' 



Figure 16: English File specifying a Szegedy quantum walk operator W(M) for an 
M with Nb = 2. Right hand side of file is not visible. 



© O b£ sze_pk.txt 



I Ph— 0- 
(1-C8-+- 



-Ry 



(1 -(0 -Ry-(Z 
Ry--(Z--(l--(0 

I Ry-.(l--(0 

0— 0— Ph I 
I Ry— (1— (0 

Ry--a~Cl~C0 
(1— (0— Ry— (Z 
(1— (0— +— Ry 



Figure 17: Picture File corresponding to the English File of Fig, 16 



D Appendix: Wocjan-Abeyesinghe Algorithm 

In this appendix, we will review the Wocjan-Abeyesinghe Algorithm for quantum 
simulated annealing, for which QuSAnn generates a quantum circuit. This appendix 
follows closely Ref. j7]. 

For any inverse temperature (3 > 0, define Vp by the quantum circuit of Fig 18 



Fig, 18 also defines the parameters a = 1, 2, . . . which we refer to as the number of 
probe bits, and c = 1, 2, . . . which we refer to as the number of phase estimation 
(PE) steps. 

Let 

Q = (e^) p ° a ° (78) 

and 

R p = vM {2Nb) ®Q]V p . (79) 
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2N B +ac 

bits 



Figure 18: Definition of operator Vg used in Ref. [7j. 



Consider an annealing schedule /3q, Pi, ■ ■ ■ , Pt r For each "time" £g{0,1,. ..,£/}, 
define a unitary matrix Up t]df recursively as follows: Let 



•)(2N B +ac) 



(80a) 



and 



(80b) 



where d G {0, 1, . . . , df — 1}. See Fig 19 for a pictorial representation of this recursion. 
We will call the parameter df = 1,2,... which is the final level of recursion, the 
Grover depth. 

If U is defined by 

^ = Uf3 tf _ 1 - df ■ ■ ■ Up 2 . df Uf 3l . df Uf3 . d} , (81) 

and if we've done things right, U should satisfy 



(82) 
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Figure 19: Pictorial representation of Grover recursion Eq.(80) 
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